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1. Introduction 

1.1. Experiment and theory 

Disorder leads to important physical effects which are of quantum mechanical 
origin. This has been revealed by Anderson ^ in the study of a disordered tight- 
binding model. This problem has attracted great attention over many decades. A 
breakthrough came with the scaling theory of localization |2j. All states in a one- 
dimensional system (1-D) are localized, whereas in 3-D a metal-insulator transition 
occurs. An analytic solution is only known for the 1-D problem 0. Although there 
was no general analytical solution available, there was consensus that in 2-D all states 
are localized. 

For quite some time after the advent of the scaling theory, many believed it to 
be essentially under control. This view is less secure nowadays, in part because recent 
experiments have challenged conventional wisdom about disordered 2-D systems. The 
2-D case still presents a problem which has become apparent by experiments 
These experiments are still being discussed controversially. Experimental reality is 
certainly more complex than a simple tight-binding model, but these results provide 
a good reason for revisiting this classic theoretical problem. 

Recently we have been able to solve the 2-D case analytically |H| ■ We have shown 
that in principle there is the possibility that the phase of delocalized states exists for 
a non- interacting electron system. For energies and disorder, where extended states 
may exist we find a coexistence of these localized and extended states. Thus the 
Anderson metal-insulator transition exists and should be regarded as a first order 
phase transition. Consequently we have returned to the old idea of Mott E| that 
the metal-insulator transition is discontinuous. This alternative idea was in its history 
completely abandoned with the advent of the scaling theory of localization. There is 
now a renaissance of it. 

This result implies the failure of the scaling theory of Anderson localization. 
Although this paper is published [H] and constitutes the basis for a new analytical 
investigation of the Anderson problem in the present paper for the case of higher 
dimensions (N-D problem) one has to accept the following: (i) the paper |S] requires 
an independent confirmation, which requires a certain time; (ii) in the history of the 
problem one has developed quite a few conceptions and this leads to a critical attitude 
towards the new theory - i.e. there is a resistance - as it departs from conventional 
wisdom. This asks for a critical evaluation of the new theory and its results. One 
should acknowledge that the problem is many-sided and quite complex. Thus logical 
errors are possible which do not lie at the surface. Perhaps a reference to B. Pascal 
is appropriate here: A truth is so delicate that any small deviation from it leads you 
to a mistake, but this mistake is also so delicate that after a small retreat you find 
yourselves in a truth again. Without an exact analytic solution any discussion cannot 
lead to firm results. The main aim of the present article is thus the generalization 
of the mathematical tools of the previous article jS] to the case of higher dimensional 
spaces and a physical interpretation of the new results. 

1.2. Structure of the present article 

The outline of the present article is a follows. In chapter 2 we give a short derivation 
of equations for N-D Anderson localization problem, and necessary summary of the 
results of the article [S]. The connection between the Anderson localization problem 
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and signal theory is discussed and the most important concept of the proposed 
methode - the filter H{z) - is defined. To understand fully these aspects of the theory 
a knowledge of the first paper (2| is recommended. The filter H{z) is generalized 
for higher dimensions D of the space. The investigation of its properties and the 
corresponding physical interpretation constitute the content of several chapters. The 
theory for high-dimensional systems (D > 4) is presented in chapter 3, whereas the 
theory for low-dimensional systems is given in chapter 4. An appendix deals with the 
mathematical conditions for the physical interpretation of the filter H{z) and a more 
detailed discussion of related aspects of the problem. 



2. New methods 



2.1. Equations for correlators 

Recently we have been able to solve the 2-D case analytically ;,8:. The tight-binding 
equation in 2-D is solved for the wave function 4'n,m and the second moments 
(correlators) (V'i™): 

^n+l.m (-^ ^n,m^^n,m V^n— l,m ^ ^ '0n,m+m'5 (1) 

m' 

where the summation over m' runs over the nearest neighbours of site (n, m) in layer 
n that is in a space of dimension p = D—1 = 1. We assume taking the limit to infinite 
size L ^ oo in p-space. The equation is solved with an initial condition 

V'0,m = 0, tpi^m = a,n- (2) 

The on-site potentials £n,m are independently and identically distributed with existing 
first two moments, {Sn.m) = and (e^ = cr^ ■ After averaging the mean squared 
amplitude (V'^ = becomes independent of m (the averaging procedure that 
justifies this statement is discussed in detail in 

In the present paper we present the analytic solution for the general case D > 2. 
The knowledge of paper |HI is prerequisite for understanding the present one as it 
contains the full formalism. The generalization of this formalism to the N-D case 
presents no problem (see below). The main eqs. Q,(j2J remain valid, only scalar 
quantities become vector quantities. 

The tight-binding equation in I + p dimension is (primed indices are summed) 

V'n+l.m ^n.mV^n,m ^n-— l,m ^ ^in,in'^n,m' 7 (3) 

^m,m' — -^^m,m' ^ ^ ^m+m'',mS (4) 
m" 

(summation over m" runs over the nearest neighbours) with initial condition V-'o,m = 
and V'l.m = am- 

One introduces the correlators (the averages are taken over disorder) 

x{n)m^i = (V'«,mV'«,i) , (5) 
v{n)m,\ = (V'«,mV'«-l,l) • (6) 
From the eq.Q one easily obtains the system of equations: 

x{n + l)mj = '5m,ic^^a;(n)m,i + xin - 1)^,1 + (7) 

Cm^m'X{n)mijiLvj - £m,m'y(»^)m',l - Aa'2/(n)l',m, 

y{n + l)m,i = -y(ri)i,m + >Cm,m'a;(n)m',i- (8) 
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ri=0 

which turn the equations into 

[z - - Cr^^rn,i)X(z)m,l - x{l)^j = Cm^m' X {z)yn' J' J (H) 

zY{z)^^i = -Y{z)i^m + /:m,m'X(z)„,.,l. (12) 

Iteration of the second equation yields 

z 1 

i^(z)m,l = -Cyn,m'X{z)m>,i 5 -X{z)m,vCv,l. (13) 

z^ ~ 1 z^ — 1 

We then obtain an equation for X(z)m,i only, which can be solved by double Fourier 
expansion. However, upon averaging over ensemble of initial conditions (as described 
in Section 3.2 of [H^): averaging over translations in p-space of boundary conditions, 
Qm) such that amCtm' = rm-m', the problem is translation- invariant in transverse 
directions. We than put: 

XizUi - / ^X(z,k)e^'^(— ). (14) 

After averaging on initial conditions the diagonal correlator becomes independent of 
m: 

3^(^)m,m = (15) 

X{zU^ ^ X{z) = j j^X{z^ k). (16) 

We obtain the final equations 
(z-1) 



(^ + 1) 



w^-£^(k)]X(z,k) = r(k)+a^X(z), (17) 



or 



5(k) 2^cos(A:j), (18) 
w — , (19) 



X(z) = ff(z)X(")(z), (20) 



X(")(z) = ii±4 



dPk r(k) 



(z-l)7 (27r)P [u;2_52(k)]^ 



(21) 



H{z) (z-l)i (27r)P[u;2-£2(k)]' ^ > 

where the X^'^^z) (or Xn^) refer to the ordered system [a = 0) and the X{z) (or a;„) 
to the disordered one (cr 7^ 0). For ea. H2U|) the inverse Z-transforni gives convolution 
property [TT]: 



^ = Y,xfhn-U (23) 

1=0 
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with 



H{z) 



Eh-n 

n=0 



(24) 



2.2. Anderson localization and signal theory 

The essential point in the analysis with respect to the localized or extended character 
of the states is to make use of signal theory from electrical engineering and switch 
from an investigation of the moments Xn to an analysis of the filter functions /i„. 

In the theory of signals jjl], Xn'' (or X^'^\z)) is the input signal, and Xn (or X{z)) 
is the output. Asymptotic behaviour of the solution is completely determined by the 
filter hn (or H{z)). The concept of the system function is a general and abstract 
description of the problem of localization. Thus the filter function has to be analysed 
to obtain general results and not the multitude of signals. This has been done in 
|H] for the 2-D case. This procedure has a certain similarity with the transition to 
an operator formalism in quantum mechanics. Particular signals x„ depend on the 
initial conditions used and do not carry much physical information because of the 
unconventional normalization For the localization problem the only property that 
matters is whether a signal belongs to the bounded or unbounded class and this can 
be derived from the filter. 

The essence of localization is contained in the filter H{z). We study the filter 
H{z) with properties described by generalized Lyapunov exponents. The filter is a 
fundamental function of the disorder a only jS| . 

A filter hn is uniquely characterized by a pole-zero diagram of its image H{z) 
which is a plot of the locations of the poles \i and zeros of H{z) in the complex- 2: 
plane. We provide just a brief summary here, for more details consult [HI EH- The 
signals Xn^ and x„ are real, therefore H{z) will have poles and zeros that are either on 
the real axis, or come in conjugate pairs. For the inverse Z-transform H{z) /i„ one 
needs to know the region of convergence (ROC). Physical considerations dictate that 
only causal filters (/i„ = for n < 0) should be considered. They have ROCs outside 
a circle that intersects the pole with max |Ai|. A causal filter is stable (bounded input 
yields a bounded output) if the unit circle \z\ = 1 is in the ROC. Note that the explicit 
calculation of hn by the inverse .Z-transform is not necessary, and it is also not feasible 
analytically due to the complexity of the function H{z). Only the type of the filter - 
stable or unstable - needs to be determined. The delocahzed states (bounded output) 
are obtained by transforming the physical solutions inside the band (bounded input) 
provided that the filter H{z) is stable. Seeking for poles is quite a simple analytical 
task which gives rather general results by elementary methods. 

As an example for the general and abstract description of the problem of 
localization as stated above let us consider the following problem. It is well-known that 
disorders broadens the band; new states outside the old band arise for \E\ > Ef, = 2D. 
Are among these new states also extended states? Numerical work for D = 3 ^1 E] 
ascertains this. This is the socalled reentrant behaviour of the mobility edge: the 
change from localized to extended states and back to localized ones upon increasing 
the disorder occurs for certain fixed energies. Because in the literature there is 
no physical explanation for this phenomenon one has simply accepted these results 
without critically examining them. 

With the help of signal theory we have found a particular transformation, which 
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gives a connection between the states in an ideal system (zero disorder) and states 
in a system with disorder. Extended states (zero disorder) as input (bounded) signal 
transform into localized states (nonzero disorder) as output (unbounded) signal if 
the filter which is responsible for this transformation is unstable. If the filter is 
stable, then extended states (bounded input signal) transform into extended states 
(bounded output signal). It is known that for zero disorder i.e. outside the band, 
\E\ > Eb, there do exist only mathematical solutions which cannot be normalized. 
These correspond to an unbounded input signal. It is impossible to find a filter 
which permits a transformation of the type unbounded input signal (mathematical 
solution) into a bounded output signal (extended states in \E\ > Ef,). The reverse 
- the transformation of an unbounded input signal (mathematical solution) into a 
unbounded output signal (localized states in \E\ > E^) - is on the other hand possible. 
I.e. the mathematical procedure developed by us generates in this case a negative 
answer to the posed question. Because this result contradicts the numerical work 
| 12l 113) . it is necessary to discuss in detail the quality of the numerical work (see 
below and Appendix B). 

We have shown in |H] that the filter H{z) is a non-analytic function of the complex 
variable z; this result remains valid also in the multi-dimensional case. The unit circle 
= 1 divides the complex plane into two analytic domains: the interior and exterior 
of the unit circle. The inverse Z-transform is quite generally defined via countour 
integrals in the complex plane 

h^^ -L f H{z)z^-. (25) 
Zm J z 

and this definition is only possible in an analytic domain. In this way in the formal 
analysis of the problem multiple solutions result. The first solution H+{z) is defined 
outside the unit circle and always exists. The filter H+{z) describes localized states 
and it is possible to connect its properties with the notion of the localization length 
|S]. The second solution H-{z) is defined inside the unit circle and does not always 
represent a solution which can be physically interpreted (this is the mathematical 
consequence that the filter be causal). The filter H_[z) describes delocalized states. 
The coexistence of the two solutions was physically interpreted in |S] as the coexistence 
of two phases - an insulating and a metallic one. Then the metal-insulator transition 
should be looked at from the basis of first-order phase transition theory. 



2.3. Conformal mapping 

The p-dimensional integral on the r.h.s. of eq. H22|l can be reduced to a one-dimensional 
integral. Consider the identity 

(y + E-=i2cos(%))dy 

^ ^ -. (26) 



The integral representation of the Dirac (5-function and the Bessel function 

-^0(2:) = 7^ / e"'=°'^('=)dfc. (27) 

will be used. Following jH], let us define a complex parameter w = u-\-iv in the upper 
half-plane, v — Im{w) > 0. Using the methods of complex variable theory we get 

dk 1 , , , ^ 

—Yd{w,E), (28) 



{2n f J w^-e^k) iw 



Exact analytic solution of the multi- dimensional Anderson localization 



7 



where 



Yd{w,E)^I [Jo{2t)f^^ cos{Et)exp{iwt)dt. (29) 



Changing the complex variable z to the parameter w corresponds to the conformal 
mapping of the inner part [\z\ < I, w — — (z^/^ + z^^/^))] or the outer part [\z\ > 1, 
w — (z^/^ + z^^/^)] of the circle onto the upper half-plane. The circle itself maps 
onto the interval [—2,2]. Note also that if H{z) — has complex conjugate poles, 
then on the upper w half-plane they differ only by the sign of u = Re(w). To avoid 
complicated notations, we seek for poles in the sector u > 0, w > and double their 
number if we find any. 

The inverse function [H] 

9 

7/7 ID I 

z= + — ±|-Vw'-4 (30) 

is double- valued. Its branch with the minus sign maps the w sector onto the inner 
part of the half-circle (|z| < 1). The second branch with the plus sign gives a mapping 
onto the half-plane with the half-circle excluded (|z| > 1). Therefore in the parametric 
w-representation 



(^ + 1) 

and 



= ±^=^ (31) 



{z - 1) 



l±^^^Y£tlB. (32) 



H±{z) 



3. High-dimensional systems 

3.1. Upper critical dimension 

It is generally assumed that 2-D systems mark the borderline between high and 
low dimension ^1]. The existence of a transition in 3-D is not questioned (high- 
dimensional systems). These assumptions, however, originate from the scaling theory 
of localization. Here the marginal dimension is Dm = 2, and a phase transition exists 
only for D > Dm- Thus perturbation theory |3] for Z? = 2 e (e ^ 1) is possible. 
The effect of statistical fluctuations cause a change of regime at I?c = 4 |1U[ [T5] ; in 
this way the upper critical dimension for localization Dc arises. For D > Dc there 
should not exist a phase transition. 

On the one hand these statements referring to higher dimensions are numerically 
nearly impossible to ascertain jj^l ■ Statistics is bad and the length of the system L 
very small (see Appendix B). Even for 13 = 3 progress towards extracting reliable 
numerical estimates of critical quantities has been remarkably difficult J7]. On the 
other hand, if the results from the scaling theory of localization for 2-D systems are 
faulty (this is what we claim), then the corresponding division of the systems into 
low and high dimensional ones is also wrong. Here one must develop an alternative 
picture. 

In the theory of critical phenomena [181 [TO] many systems belong to a class, 
where an upper critical dimension Dq has a totally different physical meaning. It 
denotes the dimension, from where on the mean field approximation is exact, or 
where with other words all critical exponents reach stationary values. I.e. in this 
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case the phase transition does exist also for D > Dq; only in the limit D — > oo 
the transition disappears, simply because the corresponding critical values have gone 
to infinity. Systems of quite different physical nature may have the same value of 
the upper critical dimension Dq. E.g. one finds Dq — 4 not only for the theory 
of magnetism [181 IT^ (in this case there exists also the marginal dimension, below 
which no phase transition is possible), but also in kinetics (cooperative phenomena 
in bimolecular processes by diffusion-controlled reactions) |2D1I^; in the latter case, 
however, there is no marginal dimension. Another example is percolation, where the 
upper dimension is Dq = 6 \22\. 

3.2. Stability and poles: solution H^{z) 

Let us start first from a purely mathematical comment: the integral Yd (w, E) is always 
finite for all w = u -\- iv in the sector u > 0, d > for high-dimensional systems with 
D > 4. This fact clearly follows from the asymptotic behaviour of the Bessel function 
for large values of its argument, Jo{2t) « cos(2i — 7r/4). We shall further on see 
that the dimensionality D = Dc = 4 is critical for localization, although this is no 
proper upper critical dimension D — Dq, which we shall determine below. Let us 
consider therefore the properties of the filter-functions in this region of the value of 
the dimension of space. 

Let us consider first the solution iJ+(z). According to 11 , the ROC of the causal 
filter is defined by the inequality \z\ > max |Ai|, where are the poles. The case when 
the system function has a pole at z = A > 1 is the simplest one for an interpretation. 
In terms of signal theory the filter H^{z) is unstable since the pole lies outside 
the unit circle \z\ = 1 in the complex x-plane. As shown in such a filter describes 
exponentially localized states (insulating phase). In order to see this one can exploit 
a basic inverse Z-transform |11J: 



Therefore the pole of Hj^{z) at z = A = exp(27) leads to exponential growth of the 
system function /i„ which in turn implies [eq. H23|) ] an exponentially increasing mean 
squared amplitude a;„. The growth exponent 7 is the so-called generalized Lyapunov 
exponent .3^ related to the localization length by ^ = 7"^. 

The value of the Lyapunov exponent 7 defines the phase. We start from the 
mathematical definition that all states with 7^0 belong to the insulating phase. The 
states with 7 = on the other hand correspond to a metal. According to this definition 
the states with non-exponential localization also belong to the metallic phase, because 
they correspond to the value 7 = 0. We can consider these states as a bad metal, in 
contrast to a good metal, where one has truly delocalized states. 

We would like to give here a summary of the results which emerge from an analysis 
of the pole diagram (for details see the Appendix A). The function Yo{uo, E) defined 
by eq. H29|l is purely imaginary for -y = 0, u > Mq, 



and the system function H+{z) itself is real. The pole must be located (if present at 
all) exactly in this region of the parameter w. It can be found as a solution of an in 
general transcendental equation 



H{z) = z/{z ~\)^hr,^ A". 



(33) 



uo = 2p+\El 



(34) 



a^^+{u) ^ 1, 



(35) 
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where 

1 /"^ 

n+{u) ^ / [M2t)f^'^ cos{Et)sm{ut)dt. (36) 

Vu - 4 Jo 

For high-dimensional systems with _D > 4 the function f2+(M) has a maximum at 
u = Up and decreases monotonically for u > uq. The system function has a pole if the 
disorder exceeds a critical value, a > ao{E), where 

ao{E) = n+{uor'^\ (37) 

Therefore, in high-dimensional systems exponential localization takes place only if the 
disorder is strong enough. 

For a < <To{E) the function H^{z) has no poles, its ROC its \z\ > 1. This 
means that the unit circle \z\ = 1 belongs to the ROC, and the filter is stable. We 
interpret solutions for this range of the disorder values, a, also as localized, however 
with non-exponential localization (which could be a power-law). Here we encounter 
limits of applicability of the method, which come into play, however, only in physically 
inaccessible systems of high dimensionality. The elementary pole search can be applied 
only for exponentially localized states, the general case requires a detailed investigation 
of the filter. 

The curve for ao{E) (Fig. la) is the well-known mobility edge. For a > (Ja{E) 
there exists an insulating phase, whereas for a < ctqIE) a metallic phase is found 
(bad metal). A coexistence of phases is here not possible. Thus the phase transition 
here has an appearance as if it were a transition of second order. This simple idea, 
however, is contradicted by the behaviour of the Lyapunov exponent: the transition 
from 7 = to 7 7^ is not continuous. One can clearly see that all these mobility 
edges for high-dimensional systems show the same qualitative behaviour. 

3.3. Stability and poles: solution H^{z) 

Now let us turn to the second branch of the solution, H^{z). In this case existence 
of the poles leads to principally different consequences. Let us assume that the 
corresponding value of the parameter w is found and the pole z = Ai is located 
inside the unit circle, |Ai| = 1/A with A > 1. Formally, however, from the definition 
w = — (z}!'^ -|- 0~^/^) the same value of w can be obtained for z = A2 = 1/Ai = A. The 
complex number A2 lies outside the region of definition of the solution, \z\ < 1. In this 
sense the pole at z = A2 is virtual. For the inverse Z-transform this fact is, however, 
irrelevant. The ROC for a causal filter is defined by the inequality \z\ > niax|Ai| or 
|z| > A > 1. Since the ROC and the region of definition of the solution \z\ < 1 do 
not intersect, a physical solution is absent. Therefore, the filter H^{z) as a physical 
solution is acceptable only if either there are no poles or they lie on the unit circle. 
The latter case is realized for D — 2 % and corresponds to so called marginal stability 
In the following, we consider the general case D > A from a unified point of view. 
A pole of the first type is related to the singularity of the function eq. (I32|) at 
w — 2 (the root, for details see the Appendix A). We define the phase and amplitude 
via the integral 129() 

Yd{w,E) ^ lD{w,E)exp{tdD{w,E)). (38) 

It is not difficult to show that this pole emerges at arbitrarily small disorder for a 
negative phase i?£)(2,£') < 0. The equation i?£)(2,£'o) = defines the boundary of 
the region \E\ > Eq, where the physical solution is absent and, therefore, any disorder 
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Figure 1. Threshold disorder values: (a) ijo{E) for the transition from the 
non-exponential to the exponential localization and (b) a-Q{E) for the transition 
from the delocalized to the (non-exponentially) localized states. The curves are 
enumerated with the values of D. 



transforms the delocalized states into localized ones. For high-dimensional systems 
the delocalized states transform into states with non-exponential localization. The 
corresponding Eq values are Eq = 3.915 (D=4), Eq = 4.365 (D=5) and Eq = 4.578 
(D=6). 

For the region \E\ < Eq there exists the physical solution with a < a'Q{E), where 
ct'q{E) is a second threshold disorder value. The behaviour of this curve ^^{E) (Fig. lb) 
is determined by resonance phenomena. The integral in (|29|l consists asymptotically 
of a power function, t^^"^)/^, and a product of trigonometric functions. E.g. the 
function cos(2t — 7r/4) comes from every Bessel function. If we represent this product 
as a sum of monochromatic waves, then we denote the existence of a wave with zero 
frequency as a resonance. For w — Q the first non-trivial resonance lies either at 
Ec = 2 (D=4,6,...) or aXEc = A (D=5,7,...). 

For \E\ < Ec a. pole of the second type appearing at higher levels of disorder 
(7 > (Tq {E) must be considered. This type of pole emerges at purely imaginary values 
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of the parameter w = iv. It corresponds to the roots of the equation 

a2f7_(t;) = 1, (39) 

where 

1 

^^iv) = —= / [Jo(2i)]^"^cos(M)cxp(-t;t)dt. (40) 

\/4 + Jo 

A physical solution is acceptable here only if the poles lie on the unit circle, = 1 
or w £ [0, 2] for our sector w. A marginally stable solution corresponds to the value 
■u; = (or u = in eq. H4()|l). The threshold disorder value is given by 

a'^{E) = ^_{Q)-^l\ (41) 

For a < cr'^{E) the function H-{z) does not possess poles, here we find the region of 
stability of the extended states. 

Spaces of dimension _D = 4, 5 possess a certain sensitivity with respect to the 
resonance phenomena mentioned above. For E^ < \E\ < Eq the line of the poles 
exhibits a deviation from purely imaginary values of the parameter w and touches the 
real axis at the point u' G [0, 2]. This point corresponds to the condition for the phase 
of the integral doW tE) = and also belongs to the unit circle (marginal stability). 
The corresponding threshold disorder value we denote again as a'f){E): 

We see that the function CTq {E) is in general singular in the energy. In going from the 
energy value Ec to the value Eq the parameter u' increases monotonically and reaches 
finally the value u' = 2. The function (t'q{E) goes continuously to zero for \E\ Eq. 

For dimensions D > 6 there are no resonance phenomena. These resonances 
are so weak that only one equation H41|l remains valid in the whole range of energies 
\E\ < Eq. I. e. although Dc = 4 is a certain critical dimension (here non-exponential 
localization arises), a qualitative agreement of all results is only obrtained for D > 6. 

In this way it emerges from our exact analytic theory that for Anderson 
localization the upper critical dimension is Dq = 6, i.e. the problem in a certain 
way shows a similarity to percolation |22| - and this may not be totally unexpected. 
However, we also note that as a rough estimate for the upper critical dimension the 
value Dc = 4 may also be accepted. It corresponds to a first and important step of the 
qualitative saturation of the results (the possibility of existence of non-exponentially 
localized states), whereas a complete saturation only obtains at D=6. 

Quite generally in the range \E\ < Eq and under the condition a < o'q{E) there 
exists a region of existence of stable delocalized states. Because in this range also non- 
localized states coexist with the other ones, this means that in this range the good 
metal (delocalized states) and the bad metal (non-exponentially localized states) can 
form a heterogeneous system, whose properties depend on the relative proportions of 
the subsystems; e.g. the bad metal, if the subsystem of non-exponentially localized 
states percolates. Fig. Q] presents results of a numerical solution of the resulting 
equations. 
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4. Low-dimensional systems 

4-1. Analytic solution 

Independently from whether Dc = 4 or Dq = 6 are taken as the upper critical 
dimension, it is necessary to consider 3-D systems as low dimensional ones. We 
approach here a widely held opinion that all states in D = 3 should be stable against 
perturbations and should have a finite radius of convergence. This, however, is not 
correct and would only be valid if the whole field of the scaling theory of localization 
would be faultless. 

The ea. (|35|) has always a solution for the physically important cases I? = 2,3. 
We refer to these cases as 'low-dimensional' ones. For the low-dimensional systems 
the integral eq. (|3(i(l can be evaluated analytically (the formulas can be found in the 
tables of Laplace transforms). The corresponding pole diagrams are shown graphically 
in the Appendix A. The system function has a pole at z = A = cxp(27) with 7 > 70, 
2 cosh(7o) = uq = 2p-\- \E\ and p = D — 1. Note, however, that the above mentioned 
feature of the low-dimensional systems is caused by the divergence of the integral 
fl+{uo). 

If one applies the obtained equations l|41|) for ctq (E) to low-dimensional systems 
with D = 2,3, then formally ao{E) = 0. Therefore, even infinitesimal disorder leads 
to solutions with exponential localization. Because the curve (To{E) forms the border 
between exponentially localized and non-exponentialy localized states, this simply 
means that in this model for low-dimensional systems non-exponentially localized 
states are impossible (see, however, a discussion below). 

For such systems the second curve (Tq{E) takes over the role of the mobility edge. 
The shape of this curve (Fig.2.b) for low-dimensional systems has a certain similarity 
with the ones for high-dimensional systems. Here we must stress that cr'^{E) does 
not represent a true mobility edge, but an upper limit for the coexistence of the two 
phases. 

For D = 2 in 8 it has been found analytically that cr^(£') = 2(1 - E'^/Af^^ for 
E < Eq where Eq = 2 = D. I.e. in this case the value Eq corresponds exactly to 
the midst of the band half-width 2D = 4. For E > Eq any disorder transforms the 
delocalized states into localized ones. 

For D = 3 the value Eq = 3 = D corresponds again to the midst of the band 
half-width 2D = 6. The resonance value mentioned above Ec lies for D = 3 at £'c = 4. 
Because Eq < Ec, the resonance does not play a role here. The mobility edge cr'aiE) 
again follows from eq. (|41|l . 

The case of the 3-D system shows an exception in the energy range Eq < \E\ < 
Eq = 3.367, where the function (Tq(£') exhibits so-called reentrant behaviom-. The 
proof of this requires a special investigation (see pole diagrams in the Appendix A). 

Note that for the 3-D system in the band center (Jq{0) = 0, because the integral 
in eq. H4()|l diverges. Fig. |21 presents results of a numerical solution of the resulting 
equations. These results are compared to the analytical result of |S] for the two- 
dimensional system. 

4-2. Phase diagram and logarithmic divergence 

Although the results for D = 3 are formally exact, they require a special discussion. 
The mathematically correct results are physically acceptable only if they are stable 
against small perturbations (e.g. small changes in the model definition, fluctuations of 
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Figure 2. Threshold disorder values: (a) ao{E) for the transition from the 
non-exponential to the exponential localization and (b) 0"g(-E) for the transition 
from the delocalized to the localized states. The curves are enumerated with the 
values of D. The curves 3' and 3" correspond to the model ea. l43l for k = 0.01 
and K = 0.1. 



parameters in the model). From this point of view two results have to be questioned: 

(i) the singularity of (t'q{E) in the band center, (7o(0) = 0, and (ii) the non-existence 
of non- localized states, <Jq(E) = 0. Point (i) arises mathematically as a consequence 
of the logarithmic divergence of the integral H4U|) for w = and \E\ 0. Point 

(ii) arises again via a logarithmic divergence of the integral l|36|l for u — uq. I.e. 
basically, the case _D = 3 is nothing else but a type of logarithmic deviation from the 
high-dimensional case. Under certain conditions one can find the perturbations which 
are capable of regularizing the mentioned logarithmic divergence, i.e. to transform 
them into a finite term. One could e.g. surmise that in this regularization correlated 
disorder |23) might play an important role. One must stress here that we are dealing 
only with results for tight-binding Hamiltonians with diagonal disorder. It is largely 
unclear, whether this property remains valid also for non-diagonal disorder. 

To illustrate this let us first consider a purely mathematical problem. Is it at all 
possible to confirm the results (i) and (ii) either by different numerical computations 
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or analytical approaches? The answer is - no! We have pointed out (Appendix B) that 
any deviation from the exact solution (via numerical or analytical ways) automatically 
generates results of a mean-field theory. A mean-field theory has a certain qualitative 
agreement with the results of the exact theory, but only for high-dimensional systems. 
Everything looks as if any uncontrolled deviation from the exact theory has added 
additional dimensions to a 3-D system. As a mathematical model let us further 
consider a 4-D system, where, however, a coupling involving this additional dimension 
is exceedingly week (parameter k <^1). We start from a generalization of the function 
(UHJ) for D = 3: 



Here the term involving n corresponds to the hopping matrix element into the fourth 
dimension. Fig[51 gives the numerical results for 2 cases: k = 0.01 (curve 3') and 
K, = 0.1 (curve 3"). One clearly sees that all results, which have nothing to do with 
the logarithmic divergence, e.g. the entire curve cfq{E) with the exception of the point 
E = 0, remain extremely stable. Even the change in the parameter Eq lie are of the 
order of O(k^). 

The 'logarithmic' results on the other hand turn out to be completely unstable. 
Even the small values of the coupling (or regularization) parameter k produce results 
which are in qualitative agreement with those for high-dimensional systems. The 
curves 3' and 3" look as if they were an extrapolation of the corresponding curves for 
ao{E) from Fig. Here the inequality ao{E) > (Jq{E) also apphes. Reducing 
the parameter k slowly moves the mobility edge oq^E) downwards, because the 
dependence on the parameter k is extremely week (logarithmic), ao{E) ~ l/ln(«;~^). 
Only for even smaller values of the parameter k <C 0.01 one might perhaps see traces 
of the exact results, because in this case one has <to{E) < (j'q{E) for \E\ < Eq. 

5. Conclusion 

Although the formal investigation of the Anderson model of localization (tight-binding 
Hamiltonian with diagonal disorder) for higher spatial dimensions D might at first look 
very abstract, the exact analytical results supply us with clear physical consequences. 
The Anderson problem of localization and the percolation problem belong to the same 
class of critical phenomena: both have the same lower and upper critical dimensions. 
I.e., although the Anderson model appears to be much more complex and richer, 
certain fundamental results appear to be transferable. Percolation is possible for 2-D 
systems, this corresponds to the existence of a metal-insulator transition in disordered 
2-D systems. In this sense there is no reason to believe that the existing contradictions 
between theory and experiment for 2-D systems point to an incompleteness of the 
Anderson model. On the contrary, our analytical investigation has shown that a 
tight-binding Hamiltonian is presumably sufficient for this purpose. 

The main problem of the theory thus does not rest in the Hamiltonian, but 
rather in the interpretation of the results, which mainly derive from numerical work. 
Our analytical and exact results demonstrate the necessity of interpreting the phase 
transition in the framework of first order phase transition theory and this holds 
independently of the spatial dimension D > 2. If one, however, attempts to apply 
procedures which have only been developed for systems with a second order phase 
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transition (and this is the general case), one does not necessarily obtain wrong 
numbers, but an incomplete or even wrong interpretation. See an example in 
where an analytical (exact) scaling function for the 2-D system has the same form 
as obtained by numerical scaling. The physical interpretation is, however totally 
different. Concerning this point we have supplied in the present article ample material 
for discussion (Appendix B) . We hope that the necessary corrections of the numerical 
tools are possible to detect the first order phase transition. 

The rather large value of the upper critical dimension for the Anderson 
localization (and the percolation) problem permits to consider 2-D and 3-D cases 
as low-dimensional systems. Thus a revisiting of the results is also necessary for 
3-D systems. We have found that the 3-D case is nothing else but a type of 
logarithmic deviation from the high-dimensional case. As a consequence results a 
certain instability of the results, whose details are discussed in the text. 

We give here a short summary of the main results. 

• For the Anderson localization problem there exists an upper critical dimension 
Dq — 6. This value is also characteristic of a related problem: percolation |22| . 
For D > Dq all phase diagrams are qualitatively the same, only the corresponding 
critical values develop in a monotonic way. One can also say that this is the 
property of a mean-field theory, although in this case a mean-field theory does 
not exist as a closed theory. 

• There exists also a second upper critical dimension Dc = 4, which has a different 
meaning. The states with non-exponential localization are formed only for 
D > Dc, whereas for D < Dc localized states are always exponentially localized. 
This second upper critical dimension Dc divides the dimensions into two classes: 
high dimensions with D > Dc and low dimensions with D < Dc- 

The results for the nontrivial spatial dimensions D > I can be summarised as 
follows. 

(i) All states with energies \E\ > Eq are localized at arbitrarily weak disorder. 
The value of Eq depends on D and lies inside the band Eq < 2D. 

(ii) For \E\ < Eg states are only localized if the disorder cr exceeds a critical 
value (To(-E), otherwise a two-phase system is formed from an insulating and a 
metallic one. This differs from the traditional point of view which considers the 
localization transition as a continuous (second order) transition. Should the standard 
interpretation of this system in the framework of first-order phase transition theory 
be applicable (which still has to be investigated) one can expect that qualitatively it 
has similar properties as other heterogeneous two-phase systems (e.g. the coexistence 
of water and ice). Then percolation problems might be important. 

(iii) Exponential localization always exists for the physically important cases 
D = 2, 3. Non-exponential localization occurs only for higher dimensions D > Dc = 4. 
In this case for < Eq and a < a'Q{E) first the heterogeneous system appears, where 
the difference between the two phases may be small (this has to be investigated). 
For ao{E) > a > ct'q^E) one finds a homogeneous system with non-exponential 
localization. Only for a > <tq{E) a system with exponential localization appears. 

(iv) ct'q^E) is in general not an analytic function of the energy E. There exist 
certain resonances. 
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Appendix A. Pole diagrams 

Parametric representation of the pole diagram 

A filter /i„ is characterized by a pole diagram of its image H{z). The principal 
definition ea. lf^ . together with the two other ones, eqs. . (|ST|) . supply us with the 
parametric w-representation of the pole diagram, eq. H32f) . Let us rewrite this relation 
into the form 

^^^l-a^R{w,E), (Al) 

where E) is generally a complex function of the complex variable w and energy 
E. The main idea is quite simple. The fmiction H{z) has its poles where 

a'^R{w,E) = l. (A.2) 

An elementary requirement for this is the condition on the argument of the complex 
function 

argi?(w,£;) = 0, (A.3) 

because a is positiv. For the given energy value E ea. (jA.3|l defines one or more lines 
(pole lines) in the complex w-plain. As we have already discussed in the text, for 
reasons of symmetry it suffices to analyse only a sector w = u + iv with u > Q and 
t; > 0. Let us analyse one of these pole lines. For each point w on this line the position 
of the pole is determined via eq.JSOI)- The corresponding value of the disorder a is 
found from ea. (|A.2|l : 

a ^ R{w,E)-^/'^. (A.4) 

Because for each value of w on the pole line there is associated a value of ct, it is 
possible to indicate by an arrow next to the line in the diagram in which direction the 
poles move with increasing disorder. 



Filter H+{z) 

Let us consider first the simplest case, the filter H{z) = H+{z). Here the poles are 
connected with the notion of the Lyapunov exponent. In the parametric representation 
there exists a simple relation 

w = 2cosh(7); (A.5) 

thus the eas. (|A.4|l . (|A.5|l together supply a connection between a and 7. 

For a 2-D system it is possible to analytically evaluate the corresponding function 
R{w,E) 0: 

R{w, E) = 1^ [^= ^ + ^ ], (A.6) 

2Vw2 - 4 ^{w + EY - 4 ^{w - EY - 4 
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or 



R{w,E) 



1 

2Vw2 - 4 
1 



(A.7) 



+ 



\/{w + E + 2){w + E-2) ^{w - E + 2){w - E 

We clearly see that the function (|A.7p possesses for the given energy several points, 
where it diverges. These values w — —2 — E , w = 2 — E , w = E — 2 and w = E + 2 are 
real. They are nothing else but the resonances discussed in the text. To detect these 
we investigate the w-parametric representation of the integral ea. H29|) . Every Bessel 
function Jo{2t) contributes asymptotically a trigonometric function, cos(2t — 7r/4). 
In addition there exists another energy dependent trigonometric function, cos{Et). 
If we represent all these functions via complex exponentials, we obtain under the 
integral in eg. 129(1 asymptotically a product of the power function and a sum (with 
well determined coefficients) of exponents exp[j(u; — 'Wj)t]. These are the resonance 
values Wj . In the 2-D case we have a strong resonance (the function R{wj , E) diverges) . 

It is easy to establish that the function R{w,E) satisfies ea. (|A.3p for the real 
value of the parameter w under the condition 

w > uq = maxjwj}, (A. 8) 

where for the 2-D system uo = 2 + \E\. I.e. here exists a line of poles which starts 
from the point uq = 2 + \E\ (it corresponds to the value a — 0, because here 
R{uo, E)~^^^ = 0) and moves on with increasing value of the parameter a always 
along the real axis. It is also easy to find that in this case no other pole lines exist. 
This result can be generalized for higher dimensions. So one can establish that the 
value of Wo found from ea. (|34|l corresponds precisely to the condition ea. (|A.8p : we 
always have to deal with the same resonance. 



(w) 



a) Low-dimensional systems 



(w) 



b) iHigh-dimensional systems 



Figure Al. Parametric representation of the pole diagram for unstable filter: a) 
low-dimensional systems with D = 2,3; b) high-dimensional systems with D > 4. 



The structure of the pole diagram is nearly identical for all dimensions D, 
Fig lAll There exists only one pole line, which originates at the point w = uq. The 
direction of increasing a is denoted by an arrow. The only difference is that for 
low-dimensional systems the starting point uq corresponds exactly to ct = (this is 
marked in the diagram with a black circle), for high-dimensional systems the value 
CTo(-£') = R{uq: E)~^/'^ is finite (the point w = uq is marked with a white circle). In the 
latter case it is impossible to find for smaller values of the disorder a < (Jo{E) a point 
in the diagram, where the condition ea. l|A.2|l is fulfilled. The filter H^{z) possesses 
no poles. The border thus found via ea.l|^ has been defined as mobility edge. 
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Filter iJ_ (z) 

This filter has a physical interpretation only under the condition |2| that the function 
H^{z) has either no poles or that they belong to the unit circle, jzj = 1. In the 
parametric w-representation the unit circle (in the sector w — u -\- iv with u > 0, 
V >Q) corresponds to an interval on the real axis u £ [0,2]. I.e. it is necessary first 
to find the single points or even lines in this interval u G [0,2], where the condition 
eg . 11 A. 3(1 holds. Later on one must also find out, in which way the pole line leaves this 
interval. 

It is easy to ascertain that one of the possible pole lines always lies on the 
imaginary axis v, and in the interval v > 0. On this line (trivial pole line) only 
the point v — could have a physical interpretation, because the point, w — 0, also 
belongs to the unit circle. The existence of other pole lines strongly depends on the 
space dimension D and the energy E. 

As the simplest example let us first consider the filter for the 2-D system |S], 
where the corresponding function has an analytical form: 

R{w, E) = / [-^ ^ + ^ =]. (A.9) 

2V4 - w2 - (w + ^A-{w-EY 

Because the interval u € [0, 2] plays an important role for the physical interpretation, 
we consider first the resonance values {wj}, which lie precisely in this region. It is easy 
to calculate that for each energy value only one resonance w — uq exists. For energies 
< < 2 we have mq — 2 — In the remaining energy region, 2 < \E\ < 4, 
one has uq — \E\ — 2. In the 2-D system the resonances are so strong (the function 
R{w,E) diverges at w = uq) that these determine in a unique manner the beginning 
of the new pole lines and their direction, Fig lA2l 




Figure A2. Parametric representation of the pole diagram for 2-D system: a) 
energy range < \E\ < 2; b) energy range 2 < \E\ < 4. 

In the interval < < 2 the pole line starts exactly at the point uq — 2 — \E\. 
This point corresponds to the value of the disorder a = and is denoted by a black 
circle in the diagram. With increasing value of a the pole line follows exactly the real 
axis u, until it reaches the point u = {w = 0). At this point the corresponding 
disorder R{0, E)^^^^ is finite (white circle in the figure), and we have denoted this 
in ea. (|41|l as aQ{E). Because this line is always a unit circle (marginal stability), we 
have physical solutions. These exist only under the condition that < a < a'Q{E). If 
the disorder crosses the border a'Q{E), the pole line leaves the point w = and follows 
further the trivial pole line (imaginary axis). Here, however, a physical interpretation 
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is no longer possible and consequently there are no extended states in the range 

In the interval 2 < \E\ < 4 the behaviour of the pole line determines the resonance 
at w = uq = \E\ — 2. This line again starts at the value uq, which corresponds to the 
disorder a = (black circle in the figure). Then this line leaves the real axis and takes 
its course into the complex plain until it reaches a particular point vo (bifurcation 
point) on the imaginary axis v. The later path makes use of branches of trivial pole 
lines. The directions belonging thereto are again marked with an arrow. Thus we have 
for cr > no points in the phase diagram which can be interpreted physically. I.e. 
an infinitesimal disorder already suffices in this energy range to destroy all extended 
states. One may also write that in this range cr'^iE) = 0. 

Formally among the 3-D case (low-dimensional system) and the cases D > A 
(high-dimensional systems) there is a quantitative difference. For high-dimensional 
systems one does not find a divergence of the function R(w, E) at the resonance 
w = uo, but only a jump in the argument of this complex function. For the 3-D 
system this divergence exists, on the other hand, but in contrast to the 2-D system 
the resonance is very weak (logarithmic divergence). In all these cases the resonances 
have no direct influence: they do not generate a pole line emerging from the point uq. 
There is, however, an indirect influence of the resonances, because every resonance 
defines the argument of the function R{w, E). 




Figure A3. Different types of pole diagrams for the systems with spatial 
dimensions D > 3. For details see the text. 

In fig lA3l different types of pole diagrams are presented. Fig lA3b is typical for 
small values of the energy. Although the resonance &t w = uq exists (this point is 
marked in the figure with a black square), the weak resonance does not generate a new 
pole line. There remains only the trivial pole line, where only the point w = has 
physical significance. The value i?(0,i?)~^/^ defines again the function ctq{E) already 
mentioned. Physical solutions exist only for a < (Jq{E) (no poles on the unit circle). 



Exact analytic solution of the multi- dimensional Anderson localization 



20 



Fig lA3b is typical for large values of the energy. In the parametric representation of the 
filter function, ea. H32|) . one can in addition see a point, where the function R{w,E) 
diverges. This is the point w = 2 (a type of energy independent resonance, which 
corresponds to the factor V4 — w^) . One can see that this point takes over the role of 
the resonance at uq. The pole line starts from the point w = 2 (this point corresponds 
to the value a = and is therefore again marked with a black circle) . The continuation 
is similar to Fig lA2b : the line continues to the imaginary axis and reaches it at the 
point V ~ vq. Although the points w = uq and w = 2 are not identical, Fig lA2b and 
Fig lA3b are qualitatively similar, and they have the same physical interpretation. In 
this range of the energy an infinitesimal disorder destroys all extended states. 

An indirect influence by the resonances consists in the fact that one can change 
in their neighbourhood the argument of the complex function E) in such a way 
that the condition ea. (|A.3|l becomes valid at a point w — u' , see FigEat and FigESh- 
The point u' is not a resonance, here the function R(u' , E) remains finite (white circle 
in the pole diagram). 

Let us consider further on two examples. In the 3-D case on can divide the 
energy values into three ranges. In the range < l-E-l < 2 the resonance is found at 
Uq = \E\, which, however, has no importance. We have the pole diagram of Fig lA3li . 
In the range 2 < < 4 the resonance occurs at uq = 4 — \E\. Here the above 
mentioned point u' arises, but under the condition \E\ > Eq ~ 3.00. I.e. in the range 
2 < \E\ < Eq the pole diagram of Fig lA3b remains valid. In the range Eq < < 4 
we find a different typ of diagram, Fig lA3h . The point u' corresponds to the condition 
< m' < Uq. For E ^ 4 the point u' moves towards the value u' — 0, and thus arises 
the bridge between the new pole line and the trivial pole line. After this the diagram 
is qualitatively the same as Fig lA3b : there are no extended states. This is also valid 
in the range 4 < < 6, where uq = \E\ — 4. 

Figure IA3b describes a complicated case which does not have an unambiguous 
interpretation. Formally this is the only diagram which has two pole lines. The 
other ones consisted always of a single pole line, although bifurcation points were also 
possible. For infinitesimally small disorder the resonance at w = 2 is important, as in 
the diagram of FigElt). Here with increasing disorder cr, the pole line which starts at 
the point w — 2 (black circle), leaves the interval u G [0,2]. The corresponding poles 
have no physical interpretation, which corresponds to the annihilation of extended 
states via infinitesimal disorder. If the disorder increases further this pole line 
approaches again the interval u £ [0,2], and reaches it at the point u' (Fig. IA3fc ^ . 
The corresponding value of the disorder cti = i?(it', E)~'^/'^ is finite (white circle). The 
physical interpretation now depends to which value of the disorder 02 — i?(0, i?)~^/^ 
corresponds the point w = (white circle) on the trivial pole line. 

If (Ti < (72, then a sort of gap arises in the disorder in such a manner that in 
the range cti < tr < (72 no values of the parameter w correspond to the pole line. I.e. 
the filter H(z) has no poles in this range. Consequently a physical interpretation of 
the solution is possible here and we obtain the reappearance of extended states at 
finite disorder values in this special interval. This condition, cti < a2, is valid for 3-D 
systems only in a very narrow range of the energy, Eq < \E\ < Eq, where Eq = 3.367. 
The behaviour of the curve 'j'q{E) is shown in FigEl (so-called reentrant behaviour). 
We can clearly see that to each energy value in the range Eq < \E\ < Eq are associated 
three values of a'Q{E); these are in particular (t'q{E) = - localization via infinitesimal 
disorder, (t'q{E) = <ti - reappearance of extended states, and (t'q{E) = 02 - again 
localization. 
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Outside this range, E'q < \E\ < 4, one has CTi > (T2- For this condition there 
are no points on the pole hne which permit a physical interpretation. I.e. after this 
annihilation of extended states via infinitesimal disorder it is not possible for extended 
states to reappear. 

In the 4-D case we find a different sequence of resonances. In the range < \E\ < 2 
the resonance occurs at uo = 2 — \E\, and the diagram Fig lA3h , is valid. In the range 
2 < |£'| < 4 the resonance is found at uq = \E\ — 2. Here the point u' {uq > u' > 2) 
arises. The pole diagram corresponds to lA8H . The pole line starts from a white circle, 
which corresponds to the value (t'q{E) > 0. Only this point in the pole diagram has 
a physical interpretation, and it means the existence of extended states under the 
condition of small values of the disorder parameter, a < ct'qIE). For £' — > 4 the values 
u' and Uq move towards w = 2. Because always u' > uq, this means that the point 
u' reaches the value w — 2 before the point uq, and this occurs for Eq = 3.915 < 4. 
For E > Eq the diagram IA3b is still valid, although the series of resonances is not 
yet exhausted: a resonance arises at wq = 6 — \E\ for 4 < jE'l < 6, and another one at 
Uo = — 6 for 6 < \E\ < 8. Here also extended states are not possible. 

Appendix B. Discussion 

Exact solution and scaling theory of localization 

The scaling theory of Anderson localization uses the conductance g as the order 
parameter. It is supposed that the system size dependence of the conductance is 
determined only by the value of the conductance itself [21 ■ This relation contains 
no information about the microscopic structure of the model, and the tight-binding 
Hamiltonian is not used in the derivation of this relation. However, one should not 
forget that the scaling theory of localization constitutes a typical phenomenological 
theory. There is no reason to believe that a phenomenological theory correctly 
reproduces all critical properties of a microscopical model (in our case the tight- 
binding Hamiltonian). Between the scaling theory and the tight-binding model there 
might exist a similar relation as e.g. between the phenomenological Landau theory 
of phase transitions and the microscopic Ising model |18[ I24| . In the best case a 
phenomenological theory is in the position to provide a qualitative description of 
the phenomena. Also qualitative differences are possible. Quantitative differences 
always exist and are unavoidable. In the case of the scaling theory of localization it 
is quite possible that the assumption that the function in the scaling relation is an 
analytic function of the conductance plays a similarly critical role as the corresponding 
assumption with respect to analytical properties of the thermodynamic potentials in 
the Landau theory. The role played by exactly solvable models in the theory of phase 
transitions [23] clearly demonstrates that exact results for microscopic models can 
never be replaced by phenomenological theories. The aim of the investigation of the 
tight-binding Hamiltonian is in no case only a confirmation of the results of scaling 
theory, but much more the search for possible deviations and problematic situations. 

Reports of metallic behaviour in dilute two-dimensional electron-hole systems 
|3] render it mandatory to reexamine |17j the basic methods which have been used 
in the past years, in conjunction with the scaling theory of localization. Exact 
diagonalization and level statistics for finite, L x L, two-dimensional systems fI^ 
and transfer-matrix approach for strips of oo x L 'itj have proven in a numerical 
way the failure of single-parameter scaling in Anderson localization, at least in the 
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statistical distribution of electronic wave-function amplitudes. 
Exact solution and numeric methods 

The Anderson localization problem constitutes a multidisciplinary problem. One is 
not dealing with the purely quantum mechanical consequences of disorder in solids. 
Formally the Schrodinger equation with random on-site potentials in the tight-binding 
representation is a stochastic algebraic equation, where physical meaning can only be 
attributed to certain average values. Averaging over random potentials forces us to 
consider statistical ensembles of macroscopically different systems. In this sense the 
problem is very similar to statistical physics, especially to the statistical physics of 
phase transitions, because there also a metal-insulator phase transition is analyzed. 
In the case of phase transitions one finds that the relevant parameters (in our case 
e.g. Lyapunov exponents or localization length) are not analytical functions of the 
disorder. These non-analytical functions originate from two steps: averaging over 
random potentials and taking the thermodynamic limit L — > oo. 

The confirmation of the validity of the scaling theory of Anderson localization 
derives mainly from numerical studies. Whether the numerical work is in the position 
to simultaneously take into account all aspects of the problem can be doubted and is in 
itself the topic for a discussion (see below) . The basic numerical methods to study the 
Anderson localization problem are in no case logically closed and consistent schemes, 
which permit a totally independent confirmation of the existence of a phase transition. 
One is dealing with a type of computer experiment, where one only tests particular 
hypotheses. If the set of such hypotheses is not complete, there is no reason to believe 
that the numerical results can be interpreted unequivocally. E.g. finite-size scaling 
theory |^ is only developed for phase transitions of second order. Finite-size scaling 
theory is by itself not in the position to determine independently the order of the phase 
transition. It only checks whether the scaling function is typical for phase transitions 
of second order. In the latter case it is possible to obtain via the scaling function 
critical values and exponents, i.e. the physical interpretation of the single-parameter 
scaling function in is correct, if the phase transition is really of second order. 

If the corresponding mathematics sees no phase transition of second order (and this 
is the case in two-dimensions), then this is no clear proof for the nonexistence of the 
phase transition. 

For two-dimensions there exists a typical single-parameter scaling function, whose 
behaviour one commonly interprets as complete localization 12, 27,. We have shown 
in [5] that for phase transitions of first order the identical single-parameter scaling 
function as in |12[ I27| has a totally different physical interpretation, it only describes 
the behaviour of the insulating phase. I.e. in this case not the numbers which 
the numerical works provides should be doubted but their physical interpretation. 
Numerical scaling is not capable of analysing a system consisting of two phases. 
Scaling theory of localization predicts instead that the metal-insulator transition is 
a continuous one in 3-D, and that all states are localized in 2-D. In this way a logical 
circle has been constructed. Scaling theory of localization claims that for the Anderson 
model only phase transitions of second order are possible or none. Finite-size scaling 
theory only checks this idea and interprets any deviation from a phase transition of 
second order as the nonexistence of the phase transition. 
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Statistic at the critical point 

Measurements of coherent transport in mesoscopic disordered systems showed large 
statistical fluctuations and the non-self-averaging nature of the transport coefficients 
such as the conductance, g. This phenomenon is referred to as universal conductance 
fluctuation [2H|- The physical quantities are broadly distributed. It is generally 
accepted that a description of this system requires distribution functions of the 
respective quantities [221 EOI- 

In the article |27] one has assumed that the full averaging {N — > oo) over different 
realizations of disorder can be replaced by one realization of disorder {N — 1). Today 
one is not quite sure about this any more and one does an averaging for N ~ 10, e.g. 

= 5 in I-G- there exists the claim that a very small number of realizations 

shows the full properties of the complete statistics. This claim can obviously not be 
proved and detailed numerical investigations 16 show clearly that the averaged values 
always exhibit a certain distribution. 

Verification of the scaling theory of localization via the conductance distribution 
requires a large statistical ensemble N . In practice, however, the quality of the 
statistics is determined by the value of the product NL^, which in numerical 
investigations is in its magnitude a constant depending on the computer. We should 
not forget in this context that we are dealing here with the statistics corresponding to 
a phase transition. It is well-known ,19) . that if one replaces in the statistics of a phase 
transition a statistical sum by the maximal term, then this procedure corresponds to 
an approximation - the mean field approximation. Because mean field theories are 
only correct for higher dimensions D > Dq, the results look qualitatively as if one 
would have artificially increased the dimension D. The replacement of the statistical 
sum by a small number of rather arbitrary terms is meaningless and produces no 
physical results. The best one can say about this numerical method in the Anderson 
localization problem is that it possibly considers some mean field model. 

Fluctuation and first order phase transition 

This discussion, however, is not fully sufficient. It does not suffice to assume that a 
probability distribution exists. It is more important to find out how this probability 
distribution arises physically. Only then can one detect, what physical importance this 
probability distribution possesses. Universal conductance fluctuation is really a finger- 
print of the critical phenomenon, but docs not fit into the generally accepted picture 
that the localization-delocalization transition is of second order. Phase transitions of 
second order do not require a description via a probability distribution, although the 
fluctuations at critical point are very strong and their role is well understood. If the 
physical origin of such fluctuations remains uncertain, then their formal description 
is also uncertain. There exists for the present case an experimental result which can 
play a decisive role. Hani et al.jBj have studied via direct electrostatic probing the 
spatial structure at the metal- insulator transition in two dimensions. They found 
a coexistence of localized and metallic regions associated with 2-D MIT. Optical 
investigations suggest that the 2-D electron system becomes strongly inhomogeneous: 
coexistence of two lines in photoluminescence spectra, one of which is caused from 
metallic regions and the other proves the existence of insulating islands in the electron 
system [ST] . 

From the point of view of the theory of second order phase transitions this 
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phenomenon is quite ununderstandable, but it has a very simple explanation in the 
framework of the theory of first order phase transitions. We are dealing here with a 
phase coexistence, where the relative proportions of the two phases varies. Then 
the universal conductance fluctuations are nothing else but a direct consequence 
of the existence of fluctuations in the heterogeneous phase. Because the transport 
properties are connected with percolation and their possibility depends strongly on 
the proportions of the two phases and their respective spatial distributions, it is clear 
that for different realizations of a static disorder potential the conductance of an 
otherwise identical mesoscopic conductor will significantly differ. The phenomenon of 
phase separation which originates from the existence of a first order phase transition in 
the 2-D electron system was discussed in It is certain that in the case of a phase 
transition of first order the description of heterogeneous fluctuations via probability 
distributions is formally possible. Whether these distributions have a clear physical 
interpretation is another question. 

This is actually a weak point in the finite-size scaling procedure |12[ I77| , because 
we have shown in |S] that in systems with a phase transition of first order mean values 
are always not self-averaging quantities. True self-averaging quantities are not only 
those which do not fluctuate within the statistical ensemble in the thermodynamic 
limit, but they should also have a physical meaning, and this is not the case for first 
order phase transitions. In this case a formal averaging over the statistical ensemble 
also takes into consideration an averaging over the phases, and the resulting averages 
have no physical meaning. 
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